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There are at least three ways in which a varying gravitational constant G could affect the 
, interpretation of the recent high-redhisft Type la supernovae results. If the local value of G at 

■ the space-time location of distant supernovae is different, it would change both the thermonuclear 

Q\ ' energy release and the time scale of the supernova outburst. In both cases the effect is related to 

a change in the Chandrasekhar mass Mch oc G~ 3 ^ 2 . Moreover the integrated variation of G with 
time would also affect cosmic evolution and therefore the luminosity distance relation. Here we 
investigate in a consistent way how these different effects of a varying G could change the current 
interpretation of the Hubble diagram of Type la supernovae. We parametrize the variation of G 
using scalar-tensor theories of gravity, such as the Jordan-Brans-Dicke theory or its extensions. It is 
remarkable that Dirac's hypothesis that G should decrease with time can qualitatively explain the 
observed Am ~ 0.2 mag decrease at z ~ 0.5 (with respect to a decelerating universe) and, at the 
same time, reduce the duration of the risetimes of distant Type la supernovae as recently reported. 
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Type la supernovae (SNela) are supposed to be one of the best examples of standard candles. This is because, 
although the nature of their progenitors and the detailed mechanism of explosion are still the subject of a strong 
debate, their observational light curves are relatively well understood and, consequently, their individual intrinsic 
• differences can be easily accounted for. Therefore, thermonuclear supernovae are well suited objects to study the 
Universe at large, especially at high redshifts (z ~ 0.5), where the rest of standard candles fail in deriving reliable 
distances, thus providing an unique tool for determining cosmological parameters or discriminating among different 
| alternative cosmological theories. 

Using the observations of 42 high redshift Type la supernovae and 18 low redshift supernovae (Riess et al. 1998; 
Perlmutter et al. 1999), both the Supernova Cosmology Project and the High-z Supernova Search Team found that 
, the peak luminosities of distant supernovae appear to be ~ 0.2 magnitude fainter than predicted by a standard 
decelerating universe (qo > 0). Based on this, the Supernova Cosmology Project derived = 0.28^0^2 at lcr, for 
a flat universe, thus forcing a non-vanishing cosmological constant. However this conclusion lies on the assumption 
that there is no mechanism likely to produce an evolution of the observed light curves over cosmological distances. 
In other words: both teams assumed that the intrinsic peak luminosity and the time scales of the light curve were 
exactly the same for both the low-z and the high-z supernovae. 

More recently Riess et al. (1999a, b) have found evidences of evolution between the samples of nearby supernovae and 
those observed at high redshifts by comparing their respective risetimes, thus casting some doubts about the derived 
cosmological parameters. In particular Riess et al. (1999a, b) find that the sample of low-z supernovae has an average 
risetime of 19.98 ± 0.15 days whereas the sample of high-z supernovae has an average risetime of 17.50 ± 0.40 days. 
The statistical likelihood that the two samples are different is high (5.8er). Riess et al. (1999b) also analyze several 
potential alternatives to produce, within a familiy of theoretical models, an evolution with the observed properties: 
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distant supernovae should be intrinsically fainter and at the same time should have smaller risetimes. All the families 
of models studied so far have the inverse trend: decreasing peak luminosities correspond to longer risetimes. 

On the other hand, and from the theoretical point of view, it is easy to show that a time variation of the gravitational 
constant, in the framework of a Scalar- Tensor cosmological theory, can reconcile the observational Hubble diagram of 
SNela with an open fl\ — universe^] The starting point is simple: assume that all thermonuclear supernovae release 
the same amount of energy (E). In a simple model of light curve (Arnett 1982) the peak luminosity is proportional 
to the mass of nickel synthetized, which in turn, to a good approximation, is a fixed fraction of the Chandrasekhar 
mass (Mnj oc Mch), which depends on the value of gravitational constant: Mch oc G -3 / 2 . Thus we have E oc G~ 3 / 2 , 
and if one assumes a slow decrease of G with time, distant supernovae should be dimmer. Moreover, the time scales 
of supernovae also depend on the Chandrasekhar mass. Let us elaborate on this last point. According to the analytic 
model of light curve of Arnett (1982), the width of the peak of the light curve of SNela is given by: 

where M c - } is the ejected mass and Mj nc is the incinerated mass. Within our current knowledge of the mechanisms 
of explosion of SNela both masses can be considered proportional to the Chandrasekhar mass, and therefore we have 
t oc M^ 2 or, equivalently, r oc G~ 3 / 4 . Since the risetime for distant supernovae is obtained from semi-empirical 
models, that is a template light curve which takes into account the decline rate and the width of the peak, one 
can then also assume this dependence on G for the risetime. This expression has the right properties since distant 
supernovae have smaller peak luminosities and, at the sime time, smaller risetimes, as required by observations. 

II. THE EFFECTS OF A VARYING G 

Despite the beauty and successes of the simplest version of General Relativity (GR) , the possibility that G could 
vary in space and/or time is well motivated. Its study can shed new light into fundamental physics and cosmology 
and it seems natural in Scalar- Tensor theories of gravity (STTs) such as Jordan-Brans-Dicke (JBD) theory or its 
extensions. 

To make quantitative predictions we will consider cosmic evolution in STTs, where G is derived from a scalar field 
(f) which is characterized by a function uj = u>((f>) determining the strength of the coupling between the scalar field and 
gravity. In the simplest JBD models, u> is just a constant and G ~ <j)~ , however if u> varies then it can increase with 
cosmic time so that ui = uj(z). The Hubble rate H in these models is given by: 

d\ 2 8irp 1 A us 4> 2 



H = {a) -^ + ^& + 3 + 6^- H ^ (2) 

this equation has to be complemented with the acceleration equations for a and 0, and with the equation of state for 
a perfect fluid: p — (7 — l)p and p + 3jHp = 0. The structure of the solutions to this set of equations is quite rich 
and depends crucially on the coupling function u>((j>) (see Barrow & Parsons 1996). Here we are only interested in the 
matter dominated regime: 7 = 1. In the weak field limit and a flat universe the exact solution is given by: 

G=i±|-r 1 = Go(l + ,) 1 /^. (3) 
6 + Zus 

In this case we also have that a = (</<o)' 2w+2 ^ < - 3w+4 ' ) - This solution for the flat universe is recovered in a general 
case in the limit t — > 00 and also arises as an exact solution of Newtonian gravity with a power law G oc t n (Barrow 
1996). For non-flat models, a(t) is not a simple power-law and the solutions get far more complicated. To illustrate 
the effects of a non-flat cosmology we will consider general solutions that can be parametrized as Eq.jBl but which 
are not simple power-laws in a(t). In this case, it is easy to check that the new Hubble law given by Eq.[2| becomes: 

H 2 (z) = Hi \ M (1 + z) 3+1 /< 1+ ^ + Q R (1 + zf + S^ A 1 (4) 



1 While we were writing this paper we became aware of a similar idea independently proposed by Amendola et al. (1999). 
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where €Im,&r and (Ia follow the usual relation: + — 1 an d are related to the familiar local ratios (z — > 0): 

n M = 87rG oP o/(3HZ), Q R = l/(a RH ) 2 and fl A = A/(3ff 2 ) by: 

fW = ^ ; Oa = ^ ; 0, = ^ (5) 



g \3 + 2ioJ g g 

1 1 to 
(1 + to) 6 (1 + toy 

Thus the GR limit is recovered as to — > oo. The luminosity distance = di,(z, VLm, ^a, w ) is obtained as usual from 
the (line-of-sight) comoving coordinate distance: r(z) = J dz' /H(z'), with the trigonometric or the hyperbolic sinus 
to account for curvature (Peebles 1993). In the limit of small z we recover the usual Hubble relation: y — H$r — 
z — (1 + <jo)z 2 /2 where a new deceleration go parameter is related to the standard one by: 

90 = — + on - r - (7) 
g 2(1 + LJ) 

One can see from this equation that even for relative small values of to the cosmological effect is small. For example 
for ~ 0.2 and Oa — 0.8 we have qo ~ —0.7 while go is around go — —0.4 for to ~ 1. Note nevertheless that this 
effect, although small, tends to decrease the acceleration and therefore it partially decreases the effect in the peak 
luminosity of SNela caused by an increasing G. In summary, Eq.|| parametrizes the change in G as a function of to 
while Eqs.[Q^| parametrize the corresponding cosmic evolution. 

As mentioned in the introduction, we are assuming that thermonuclear supernovae release a similar amount of 
energy E oc G~ 3 / 2 . Thus using Eq.||, we have: 

E , G n-3/2 15 , / G s 15 , / \ 

E-r^) ; M - M ° = T log ^)-i(iT^ log ( 1+z )' (8) 

were M is the absolute magnitude and the subscript denotes the local value. Therefore we have the following Hubble 
relation: 

15 

m(z)=M + 5 \ogd L {z,n M ^K,io) + — — — - log (l + z) (9) 

4(1 + to) 

which reduces to the standard relation as to — > oo. ^,From the last term alone we can see that to ~ 5 can reduce 
the apparent luminosity by Am ~ 0.2, which is roughly what is needed to explain the SNela results without a 
cosmological constant. For illustrative purposes figure [j] shows the above relation for two representative cosmological 
models, including the effects of w in d^, for to — ±5 (dotted lines) and the standard (to — oo) case (solid line). 

The effect of a varying G on the time scales of SNela can be obtained from Eq.[fL|. Since r oc G~ 3 / 4 , the ratio of 
the local time scale, to, to the faraway one is: 

/ G \ -3/4 3 

( — ) =(l + z)- — . (10) 
to' \ Go I 

and, to make some quantitative estimates, we can use the mean evolution found by Riess et al. (1999a, b). From their 
figure 1 we obtain the following widths of the light curve when the supernova is 2.5 magnitudes fainter than the peak 
luminosity: r = 45.0 ± 0.15 (at z ~ 0) and r = 43.8 ± 0.40 (at z ~ 0.5), were the errors in the widths have been 
ascribed solely to the errors in the risetimes. Thus, from Eq.|lCj we obtain to ~ 10.25^^^ (2 a errors). Therefore, a 
very small variation of the gravitational constant can account for the reported differences in the SNela time scales. 
However these limits on to should be considered as weak, in the sense that since most SNela are discovered close to its 
peak luminosity the width of the light curve is poorly determined. These values are shown as horizontal dashed (la) 
and continuous (2a) lines in Fig. Q where the confidence contours (at the 99%, 90%, 68% — solid lines — 5% and 
1% confidence level — dotted lines) in the (lo,D,\) plane for a flat fifl = universe (left panel) and in the (w,f2jvf) 
plane for the case = (right panel) are shown. 

III. DISCUSSION AND CONCLUSIONS 

In astrophysics and cosmology the laws of physics (and in particular the simplest version of GR) are extrapolated 
outside its observational range of validity. It is therefore important to test for deviations of these laws at increasing 
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cosmological scales and times (redshifts) . SNela provide us with a new tool to test how the laws of gravity and 
cosmology were in farway galaxies (z ~ 0.5). In particular, current limits on the (parametrized) Post Newtonian 
formalism mostly restrict to our very local Universe (see Will 1993). The observational limits on G/G come from quite 
different times and scales (see Barrow & Parsons 1996 for a review), but mostly in the local and nearby enviroments 
at z ~ (solar system, binary pulsars, white dwarf cooling, neutron stars) typical bounds give G/G < 1CT 11 — 10~ 12 
yr~ x , or ui > 10 — 100. However, STTs predict u> = w{4>). That is, u is not required to be a constant, so that u> 
can increase with cosmic time, u> = w(z), in such a way that it could approach the GR predictions (u> — > oo) at 
present time and still give significant deviations at earlier cosmological times. In this sense bounds from primordial 
nucleosynthesis could provide an important test. Current bounds on uj from nucleosynthesis are comparable to the 
local values but these bounds are model dependent and also involve very large extrapolations. 

Our analysis indicates that if we adopt the constraints derived from the width of the light curves of SNela then 
our best fit to the data requires ui ~ 10 (or equivalently G/G ~ 10~ n yr -1 or ~ 10% in G). This value is slightly 
smaller than some of the the current constraints at z ~ 0, but it corresponds to higher redshifts z ~ 0.5 and could be 
accomodated in STTs with ui = cu((f)) = u>(z). If this is the case, at the 2a confidence level we obtain 0.0 < fl\ < 1.0 
and the Hubble diagram of SNela poorly constrains Om < 1. At the la confidence level we obtain 0.2 < fl\ < 0.8 
and f^M J$ 0.7. If we do not take into account the restrictions derived from the width of the light curves then 
our conclusions are much weaker: the observational data and the theory can be reconciled in the framework of a 
cosmological theory with a varying G with no cosmological constant (J7a = 0) only if uj > 1.5. If we further require a 
flat J1r = universe then 1.5 < uj < 3.0 is needed. 

Obviously more work is needed both regarding other observational consequences of STTs and on the physics of 
supernovae. In particular, an improvement of our knowledge of the physics of thermonuclear supernovae would 
provide us with an unique tool to test fundamental laws of physics over cosmological distances. In addition it should 
be stressed that new observations of distant supernovae, or other standard candles, at higher redshifts [z > 1) could 
constrain even more the current limits on the variation of the fundamental constants. 
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FIG. 1. Hubble diagram for the high-redshift SNe. 
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FIG. 2. Confidence contours in the plane (ui, Qa) for a flat case £Ir = (left panel) and in the plane (ui,Qm) for the case 
S7a = (right panel). 
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